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ABSTRACT 

Near-infrared photometric observations of the Type Iln SN 2005ip in NGC 2906 reveal large fluxes 
(>1.3 mJy) in the -ftTs-band over more than 900 days. While warm dust can explain the late-time 
iiTs-band emission of SN 2005ip, the nature of the dust heating source is ambiguous. Shock heating 
of pre-existing dust by post-shocked gas is unlikely because the forward shock is moving too slowly 
to have traversed the expected dust-free cavity by the time observations first reveal the Ks emission. 
While an infrared light echo model correctly predicts a near-infrared luminosity plateau, heating 
dust to the observed temperatures of ~1400-1600 K at a relatively large distance from the supernova 
(> 10^^ cm) requires an extraordinarily high early supernova luminosity (~ 1 x 10^^ L©)- The evidence 
instead favors condensing dust in the cool, dense shell between the forward and reverse shocks. Both 
the initial dust temperature and the evolutionary trend towards lower temperatures are consistent 
with this scenario. We infer that radiation from the circumstellar interaction heats the dust. While 
this paper includes no spectroscopic confirmation, the photometry is comparable to other SNe that 
do show spectroscopic evidence for dust formation. Observations of dust formation in SNe are sparse, 
so these results provide a rare opportunity to consider SNe Type Iln as dust sources. 
Subject headings: circumstellar matter — supernovae: general — supernovae:individual: SN 2005ip 
— dust, extinction — infrared: stars 
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1. INTRODUCTION 

For nearly 40 years, core-collapse SN events have 
been c onsidered as possible s ources of dust in the uni- 
verse (jCernuschi et al.l Il967t iHovle fc Wickramasinghd 
Il970( ). More recent studies have proposed that the 
core-collapse supernovae may be the pr imary sources 
of dust in t he early universe (iTodini fc Fcrrara 2001; 

l2007fl . Several models 



iNozawa et al. 2003; D wek et all 



Todini fc Ferrara.2001lTNozawa et aUl2003l |2008[ ) suc- 
ceed in producing the large amounts of dust observed at 
high redshifts, and recent data present the first evidence 
for a supernova origin for dust in an object at z > 6 
(jMaiolino et al.ll2004f) . 

Nonetheless, direct observational evidence for dust for- 
mation in st ipernovae remains sparse, even in the lo- 
cal universe (|Meikle et al.ll2007l and references therein). 
Any newly formed dust would produce a late-time near- 
infrared ex cess in compar ison to the blackbody optical 
spectrum. iMerrill (|1980l ) was the first to detect such 
an infrared excess from a supernova (SN 1979C). Since 
then, only a handful of near-infrared excesses associ- 
ated with core-collapse events have been observed (e.g. 
Dwek fc Arendt 1992; Gcrardv et al. 2002; Pozzo et al] 



20041: lMeikl~t al.. .2006. : .Smith et al.. .2008. . and refer- 



ences therein). 

Late-time near-infrared emission, however, is not 
unique to dust formation. Several possible mechanisms 
may give rise to such an excess due to the complex nature 
of the surrounding environment of most core-collapse su- 
pernovae (see Figure [T]). (1) As the ejecta expand and 
cool, dust condenses and is radiatively heated. (2) The 
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hot gas behind the forward shock coUisionally heats pre- 
existing circumstellar gas and dust. (3) The SN peak 
luminosity heats pre-existing dust and produces an "IR 
echo." (4) Dust condenses in either the forward or re- 
verse shocks and is heated by radiative shock emission. 

The first scenario requires dust to condense from the 
ejecta as it expands and cools. The physical conditions in 
an expanding SN support this possibility. In particular, 
large abundances of the necessary elements, cooling of 
the ejecta via expansion, and dynamical ins t abilities can 
result in dust formation (iLucv et al.lll991l : iDwek etldl 
ll992tlMeikle et al.lll993l : lRoche et al.lll993[ ). In this sce- 
nario the primary energy source is the radioactive power 
released in the expanding gas, which sets the luminos- 
ity of the dust. It is possible, however, for other energy 
sources to exist, such as circumstellar interactions. 

The second scenario describes pre-existing circumstel- 
lar dust heated directly by shock interactions. The sce- 
nario is more likely to occur at later times. Assum- 
ing a s pherica l ly sy mmetric circumstellar distribution 
of dust, iDwekl ()1983f ) shows that a SN peak-luminosity 
of ~ 10^° Lq creates a dust free cavity, via photo- 
evaporation, with a radius of ~ 3 x 10^^ cm (0.1 light 
years). For typical expansion velocities (~ 5000 km s~^), 
this cavity delays the onset of emission by a couple of 
years. In addition, dust sputtering may limit the grain 
radiation efficiency so that only a small po rtion of the to - 
tal shock power is converted into radiation ()Drainelll981h . 

In the third scenario, the energy from the SN peak- 
luminosity once again creates a dust free cavity. The 
remaining shell of dust warms to high temperatures but 
does not vaporize. Due to light travel time effects, the 
thermal radiation from the dust grains reaches the ob- 
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Fig . 1. — An illustration of the environment surrounding core-collapse supernovae similar to SN 2005ip (adapted from ISmith et all 
II2OO8I )). Dust may form in a number of locations. Following the initial supernova explosion, the SN ejecta expand and cool sufficiently for 
dust formation to occur. A shocked region forms in between a forward shock in the circumstellar gas and a reverse shock in the supernova 
ejecta. If the gas is able to cool efficiently, dust condensation is likely to occur in the cool, dense shell behind the shock. Radiative cooling 
is more likely at the reverse shock because of the higher density, lower shock velocity, and the possibility of heavy element enrichment. 
Cooling at the forward shock is less likely, but the presence of clumps in the circumstellar medium, as indicated by the relatively narrow 
lines observed in SNe Iln, may allow radiative shocks. 



server over an extended period of time (jDwekl 119831 ). 
The infrared luminosity plateau occurs on year long time 
scales, which are set by the light travel time across the 
inner edge of the cavity. An IR echo's spectrum and 
temporal evolution are useful for estimating the mass- 
loss history of the progenitor star, and to investigate the 
'eometry an d composition o f the circu mstellar structure 
Dweklll983l ITqSS: Emmcring fc Chevalier 1988) . 

In the final scenario, the ejecta collide with the pre- 
existing circumstellar medium, which creates a forward 
shock. The circumstellar interaction decelerates the blast 
wave, thereby simultaneously creating a reverse shock. 
The deceleration converts some of the kinetic energy of 
the shock into X-rays and visible light. Dust grains are 
likely to condense in the cool, dense shell that forms 
behind the rad iative shocks as they undergo a ther- 
mal instability (|Pozzo et al.l 12004 ). Typically, the post- 
shock conditions surrounding reverse shocks are more 
conducive to dust condensation as the chemically rich, 
dense ejecta pass through the shock front and cool. 

The decelerated blast wave associated with the reverse 
shock is more likely to cool in sufficiently dense circum- 
stellar environments. Type Iln events typically have the 
densest environments of any SNe, as they are defined 
by the "narrow" H lines that originate from clumps in 



the pre-existing dense H-r ich shell that is moving at 
a relatively slow velocity (Schl egeii ll990( ). Therefore, 
one might suppose dust formation in the cool, dense 
shell behind the reverse shock to be associated with 
Type Iln events. Unfortunately, Type Iln events are 
rare, consisting of on ly ~2-3% of all core-collapse SNe 
(jGal-Yam et al.ll2007D . and few well-studied events exist. 
This paper presents an analysis of the observed late- 
time near-infrared emission in the Type Iln SN 2005ip 
from days 7-935 post-discovery. In fj2l we present the 
observations, data reduction techniques, and photome- 
try. In fJSj we discuss the different possible dust heating 
mechanisms and compare the SN to SNe with similar 
characteristics. The observations suggest that either an 
IR echo or dust formation in the post-shocked gas ex- 
plain the long lived IR excess, but that dust formation 
is more likely. 21 summarizes these results. 

2. OBSERVATIONS 

2.1. SN 2005ip 

Supernova 2005ip was discovered in NGC 2906 on UT 
2005 November 5 (Boles et al. 2QQB)- NGC 2906 lies at 
a distance of 30 Mpc, given a host galaxy re cession ve- 
locity of 2140 km s~^ (lde Vaucouleurs et al.lll99li) and 
Ho ^ 72 km s^^ Mpc^ (jFreedman et al.ll2001h . Earlv 
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optical spectra suggested the discovery occu rred a few 
weeks following the explosion (jModjaz et al.1 12005). The 
presence of H in the spectru m led to the classifi cation 
of the SN as a Type II event (jModiaz et al.ll2005f) . Late 
spectra of SN 2005ip showed the development of a nar- 
row H g line that is characteris tic of Type Iln super- 
novae (jSmith et al.ll2008. in prepD . Various observers ac- 
cumulated optical photometry of SN 2005ip over the first 
100 days (Teamo 2007-'^). The supernova was brightest, 
V ~ 14.8, at the time of the first observation. The fact 
that all data were obtained post-peak does not allow for 
the explosion date to be constrained. An unfiltered R- 
band light curve beginning on d ay ~10 post-discovery 
from ISmith et al.l (|2008. in prenT ) shows a linear magni- 
tude decline through day ~150 post-discovery followed 
by luminosity plateau at i? ^ 17.8 that begins no later 
than day ~300 and persists throughout the duration of 
our observations. Observations from Swift on day 461 
reveal a magnitude V ~ 18.5 and a 0. 2 — 10 keV X-ray 
luminosity of 1.6 ± 0.3 x 10*° erg s~^ (jlmmler fc PoolevI 
[2007.) . 

2.2. Near-Infrared Photometry 

Our observations of SN 2005ip consist of J-,H-, and 
Ks-hand images spanning from 2005 December through 
2008 June, days 7-935 following discovery. All observa- 
tions were made with FanCam, a 1024 x 1024 HAWAII- 
I HgCdTe imaging system on the University of Vir- 
ginia's 31-inch telesco pe at Fan Mountain, just outsid e 
of Charlottesville, VA (jKanneganti et al]|2008. in prepi ). 
Each epoch consists of fifteen minutes of integration in 
JHKs bands, which have detection limits at the lOa level 
of 0.066, 0.098, and 0.156 mJy (or 18.5, 17.5, and 16.5 
mag), respectively. Individual exposures are sky back- 
ground limited and have an integration time of either 
30 or 60 seconds. Flat-field frames are composed of 
dusk and dawn sky observations. We employed stan- 
dard near-infrared data reduction techniques in IRAF. 
The relatively small galaxy size made it possible to fit it 
into a single array quadrant. Empty quadrants were ef- 
ficiently utilized as sky exposures. Data were taken with 
the galaxy placed in each quadrant and each quadrant 
was reduced separately. Ultimately, all reduced quad- 
rants were coadded. Figure [2^ shows a typical reduced 
J-band image. 

Template subtraction minimizes photometric confu- 
sion from the underlying galaxy. While there are no 
existing FanCam observations of the galaxy prior to SN 
2005ip, IRAF's DAOPHOT PSF and SUBSTAR pack- 
ages provide the means to create and subtract an average 
PSF at the SN location for each epoch (Figure [2)d) . A 
coaddition of each epoch yields a master template that 
is deeper and smoother than any individual epoch. The 
subtraction of the template from each epoch leaves only 
the SN (Figure [2}:) . To ensure a clean subtraction, indi- 
vidual and template PSFs are scaled to match in shape 
and intensity using IRAF's PSFMATCH package. We 
performed photometry with IRAF's APPHOT package, 
implementing an 8-pixel radius aperture and a 3-pixel 
wide annulus with an 8-pixel inner radius. For magni- 
tude calibration, the SN is compared to three 2MASS 
reference stars located in the eastern part of Figure [2l 




Fig. 2. — Figure [2^ (top) shows a J-band image of SN 2005ip 
within NGC 2906 at 7 days post-discovery. On the left hand side 
are three 2MASS catalogued stars, which serve as calibration ref- 
erences. Figure [2J3 (middle) shows the artificial template produced 
using PSF subtraction, and Figure [2j; (bottom) shows the original 
image following template subtraction. An imperfect subtraction 
leaves a slight residual, which is enhanced in Figure [2]; by the 
contrast settings. Removing the underlying galaxy via template 
subtraction minimizes photometric error. 



http://www.astrosurf.com/snweb2/2005/05ip/05ipMeas.htm 



TABLE 1 

Near-Infrared Photometry of SN 2005ip 



JD 

(2450000+) 



Epoch 

(days) 



J 

(mag) 



H 

(mag) 



(mag) 



J 

(mjy) 



H 

(mJy) 



(mJy) 



3686.5 


7.0 


14.42 


0.01) 


- 


13.87(0.01) 


2.84(0.03) 


- 


1.75(0.02) 


3688.0 


8.0 


14.41 


0.01) 


- 


13.93(0.02) 


2.88(0.03) 


- 


1.66(0.03) 


3743.0 


63.0 


15.42 


0.03) 


14.91(0.05) 


14.12(0.03) 


1.14(0.03) 


1.07(0.05) 


1.39(0.03) 


3763.1 


83.0 


15.51 


0.03) 


14.71(0.03) 


13.82(0.02) 


1.04(0.03) 


1.28(0.04) 


1.83(0.03) 


3773.0 


93.0 


15.59 


0.03) 


14.66(0.03) 


13.75(0.01) 


0.97(0.02) 


1.34(0.04) 


1.96(0.03) 


3799.5 


120.0 


15.84 


0.03) 


14.66(0.03) 


13.67(0.01) 


0.77(0.02) 


1.34(0.04) 


2.11(0.03) 


3830.7 


151.0 


16.08 


0.04) 


14.77(0.03) 


13.72(0.01) 


0.62(0.02) 


1.21(0.04) 


2.01(0.03) 


3850.5 


171.0 


16.21 


0.05) 


14.81(0.04) 


13.75(0.02) 


0.54(0.02) 


1.16(0.04) 


1.97(0.04) 


3865.5 


186.0 


15.96 


0.04) 


14.75(0.04) 


13.67(0.01) 


0.69(0.02) 


1.23(0.04) 


2.12(0.03) 


3879.5 


200.0 


16.05 


0.05) 


14.78(0.04) 


13.67(0.02) 


0.63(0.03) 


1.20(0.04) 


2.10(0.03) 


4030.9 


351.0 


16.25 


0.05) 


- 


13.75(0.01) 


0.53(0.02) 


- 


1.96(0.03) 


4112.0 


432.0 


16.49 


0.06) 


15.20(0.05) 


13.87(0.01) 


0.42(0.02) 


0.82(0.04) 


1.75(0.02) 


4159.8 


480.0 


16.60 


0.07) 


15.15(0.05) 


13.89(0.02) 


0.38(0.02) 


0.85(0.04) 


1.73(0.03) 


4227.6 


548.0 


16.66 


0.08) 


15.23(0.06) 


13.96(0.02) 


0.36(0.03) 


0.79(0.04) 


1.62(0.03) 


4367.9 


688.0 


16.76 


0.08) 


15.48(0.07) 


13.92(0.02) 


0.33(0.02) 


0.63(0.04) 


1.68(0.03) 


4406.8 


722.0 


16.85 


0.09) 


15.52(0.07) 


14.03(0.02) 


0.30(0.03) 


0.61(0.04) 


1.52(0.03) 


4452.9 


767.0 


16.94 


0.11) 


15.67(0.08) 


14.07(0.02) 


0.28(0.03) 


0.53(0.04) 


1.46(0.02) 


4516.8 


831.0 






15.87(0.11) 


14.14(0.03) 


- 


0.44(0.05) 


1.37(0.04) 


4537.6 


852.0 


17.03 


0.10) 


15.80(0.09) 


14.10(0.02) 


0.26(0.02) 


0.47(0.04) 


1.42(0.03) 


4580.6 


895.0 


17.01 


0.10) 


15.83(0.09) 


14.12(0.02) 


0.26(0.02) 


0.45(0.04) 


1.40(0.02) 


4620.6 


935.0 


17.17 


0.12) 


16.14(0.15) 


14.18(0.05) 


0.23(0.03) 


0.34(0.05) 


1.32(0.06) 
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Fig. 3. — JHKs light curves for SN 2005ip showing magnitude (top) and flux (bottom). The ii's-band curve remains fairly constant 
until late-times, even as the J-band declines. The observed kink in the H and K light curves at day ~50 is likely due to the onset of dust 
formation in the post-shocked gas. 
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TABLE 2 
Colors of SN 2005ip 



Epoch 


J-H 


J-Ks 


H-Ks 


(days) 


(mag) 


(mag) 


(mag) 


7 


- 


0.55(0.02) 


- 


8 


- 


0.48(0.02) 


- 


63 


0.51(0.06) 


1.30(0.04) 


0.78(0.06) 


83 


0.81(0.04) 


1.69(0.03) 


0.89(0.04) 


93 


0.93(0.04) 


1.85(0.03) 


0.92(0.03) 


120 


1.18(0.05) 


2.17(0.04) 


0.99(0.03) 


151 


1.31(0.06) 


2.36(0.05) 


1.05(0.04) 


171 


1.40(0.06) 


2.47(0.05) 


1.07(0.04) 


186 


1.21(0.05) 


2.30(0.04) 


1.09(0.04) 


200 


1.27(0.06) 


2.38(0.05) 


1.11(0.04) 


351 


- 


2.50(0.05) 


- 


432 


1.29(0.08) 


2.62(0.06) 


1.33(0.05) 


480 


1.45(0.08) 


2.72(0.07) 


1.27(0.05) 


548 


1.43(0.09) 


2.70(0.08) 


1.27(0.06) 


688 


1.29(0.11) 


2.85(0.08) 


1.56(0.07) 


722 


1.35(0.11) 


2.84(0.09) 


1.49(0.07) 


767 


1.29(0.13) 


2.89(0.11) 


1.60(0.08) 


831 


-15.87(0.11) 


-14.14(0.03) 


1.73(0.12) 


852 


1.24(0.14) 


2.94(0.10) 


1.70(0.09) 


895 


1.18(0.13) 


2.90(0.10) 


1.72(0.09) 


935 


1.04(0.19) 


2.99(0.13) 


1.95(0.16) 



This template subtraction method yields results sim- 
ilar to, but less noisy than, other techniques, such as 
aperture photometry and PSF photometry without tem- 
plate subtraction. Table [T] lists the JHKg photometry 
and Figure [3] plots the light curve. The i^T^-band magni- 
tude remains fairly constant at < 14.2 mag (a flux of > 
1.3 mJy) for the duration of the observations, even as the 
J-band declines. A notable 'kink' occurs at day ~200, 
although we do not speculate on the cause in this paper. 
Table [5] lists the near-infrared color evolution. Figure 2] 
shows the color to quickly redden between ~50 and 200 
days. This trend becomes more gradual and continues 
throughout the last observation epoch. 

3. ANALYSIS AND DISCUSSION 

3.1. Temperature and Luminosity Evolution 

Figure [5] shows the near-infrared fluxes. Matching the 
H — Kg colors to the Planck function, presuming a dust 
emission/absorption efficiency, (5(A), establishes a near- 
infrared dust temperature. The flux of a single dust par- 
ticle of radius a at a temperature T^ and wavelength A 
is given by 

FA,T-Q(A)i?A(T), (1) 

where B\{T) is the Planck function. The emis- 
sion/absorption efficiency is given as 



for a ^ A 



(5 = 1, for a > A 



(2) 
(3) 



where n is an integer. For large grain sizes (e.g. a > 
1 ^m), n = because the dust acts as a blackbody. For 
smaller grains (i.e. a ^ 1/^ni), n > 0. We consider the 
cases n=0, 1, and 2 to account for bot h large grains that 
are associated with condensing dust (jPozzo et al.ll2004 ) 
and smaller grains that are associated with pre-existing, 
optically thin dust, which is charac terized by A~^ an d 
A^^ emission/absorption efficiencies (|Dwek et al.lll983l ). 
Figure [5] plots the corresponding Planck curves over 
the original fluxes. The fits exclude the J-band fluxes be- 
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Fig. 4.— The J-Ks and H-Ks color evolution for SN 2005ip. 
The J — Ks color reveals the quickly declining photospheric com- 
ponent in the J-band at early times. The colors continue to redden 
throughout the observations as the warm dust dominates the in- 
frared emission. 



cause of SN photospheric flux and line emission contam- 
ination, especially at early times when the photospheric 
temperature is ~10,000 K. At these early times, the H- 
band flux may also contain some contribution from the 
SN photosphere. By later times, however, we presume 
the photosphere has faded sufficiently so that the warm 
dust dominates the flux. With only 2 data points for 
each epoch there is no clear 'best' fit, so it is necessary 
to consider all cases. 

Overall, the Planck peaks trend towards longer wave- 
lengths/cooler temperatures. Table [3] lists the H — Ks 
color temperatures for each emission/absorption effi- 
ciency at each epoch. Figure [6^ plots the color temper- 
ature evolution. In all cases, the fading SN photosphere 
contributes to the dramatic temperature decline over the 
first ~150 days. Although the exact point at which 
the dust temperature begins to dominate is unknown, 
we presume the photospheric component to decline suf- 
ficiently below the dust temperature after ~150 days, 
at which point the temperature begins a slow, steady 
decline. While the chosen efficiency can alter the de- 
rived dust temperatures at early times by a few hundred 
Kelvin, these temperatures are less than or equal to the 
dust vaporization temperature (~ 1400 — 1600 K). Over 
the next 2 years, the temperatures drop to ~1000 K. 
Considering the observations are made at near-infrared 
wavelengths, they probe the warmest and most luminous 
grains. Cooler dust may also be present, but its emission 
peaks at wavelengths longer than those observed here. 

In contrast to the decreasing dust temperature, the 
dust luminosity plateaus early on and does not decline 
throughout the extent of the observations. Assuming a 
distance of 30 Mpc, integrating the Planck curves over all 
wavelengths yields a lower limit to the bolometric dust 
luminosity, shown in Figure ^p and listed in Table [31 
Again, this calculation does not account for cooler dust 
that could contribute at longer wavelengths. Through- 
out the entire set of observations, the dust luminosity 
remains at ^ 3 x 10"*^ erg s~^. 

3.2. Comparison to Previous SNe 

Comparing SN 2005ip photometry to prior, well- 
studied events can help expose the dust heating 
mechanism. We find three well-sampled SNe with 



Epoch 

(d) 



TABLE 3 
Temperatures and Luminosities of SN 2005ip 



Ttft 
(K) 



(K) 



(K) (lO-^ierg s-l) (lO^^erg s^^) (lO^^crg s^i) 



J-'radioactive 



7 


- 


- 


- 


- 


- 


- 


42.09 


8 


- 


- 


- 


- 


- 


- 


42.09 


63 


1969(89) 


1547(56) 


1279(39) 


41.43(0.03) 


41.38(0.03) 


41.34(0.03) 


41.86 


83 


1809(47) 
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similar br ight near-infra red late-time em ission: SN 
1998S dFas sia ct al. 2000; iPozzo et al.l[200l. SN 1 995N 
(iGerardv et al.ll2002|). and SN 200 6 ic (ISakon et al.ll2007l: 



!002r 

m 



Smith et al.ri2008t iMattila et all l2008t iNozawa et all 



20081: iDi Carlo et al.1 I2008D . The associated spectra of 



these SNe reveal further details about the surrounding 
environment and the nature of the infrared emission. 

Type Iln SN 1998S is most comparable to SN 2005ip, 
with similar luminosities and temperatures, although 
a decline in the J^-band luminos i ty is observed as 
early as day 333 (iFassia et all 120001: iPoolev et"an i2002t 



early as day 666 iirassia et ai.M/uuut \foc 
IPozzo et all |2004|) . iPozzo et all (|2004[ ) 



conclude dust 
condensation in the cool, dense shell behind the reverse 
shock is responsible for the late-time luminosity. They do 
not consider the forward shock given the increased like- 
lihood for dust condensation in the post-shocked ejecta. 
All observations are consistent with this scenario. The 
observed optical line profiles likely arise from the low- 
ionization gas in the cooling shock ( Cheva lier fc FranssonI 
|1994[ ). The coincidence of the line velocities with the ve- 
locity of the dust shell (derived from the estimated shell 
radius from a blackbody approximation versus the time 
since the explosion) suggests a physical relationship be- 
tween the two. The condensing dust also explains the 
fading central and redshifted Ha peaks with respect t o 
the blueshifted pea k observed bv IGerardv et alj (|2002f ). 
iFassia et al.l (|200lD observe CO emission to develop at 
the same time the absorption is observed. Finally, the 
pure blackbody fits of the dust (n=0) are more sugges- 
tive of dust condensing in optically thick lumps than pre- 
existing, optically thin dust. 

Type Iln SN 1995N also shows com parable temper- 
atures and lum inosities to SN 2005ip JFcix et al.l 120001 : 
IGerardv et"aI1 2002; Pa storello et al.ll2005D . Similar to 
SN 2005ip, the luminosity plateau extends all the way 
through day ~ 10 00, after which a dec lin e was observed 
out to day 2413 (|Gerardv et al.ll2002D . IGerardv et all 
(|2002f ) do not consider dust condensation, but instead 
suggest the dust is heated by either an IR echo or by 
the X-rays produced in direct circumstellar interaction 



heating of pre-existing dust. While the IR echo model 
adequately describes the photometric evolution, the im- 
plied dust cavity radius, created via heating and vapor- 
ization by the peak SN luminosity, requires an unusually 
bright peak luminosity. Instead, the authors suggest an 
event prior to the core collapse could have cleared out 
the large dust cavity required in this instance. Alter- 
natively, the authors' circumstellar medium interaction 
model does not require such a large dust cavity. The 
pre-existing dust is heated by X-ray and UV emissions 
produced by the much more slowly evolving shock inter- 
actions with a dense circumstellar m edium. 

Interestingly, IGerardv et al.l (|2002f ) rule out direct col- 
lisional heating of pre-existing dust by the hot, post- 
shocked electrons. The authors argue dust temperatures 
of 800 — 1000 K require prohibitively high gas densi- 
ties {ug > 10^^ cm^^ for Tg = 10^ K). However, we 
find these densities are not required in this scenario. 
For a = 0. 1 ^tm silicate grain s and a gas temperature 
Te = lO'' K. lDwek etaLJ (|2008D show in their figure 1 that 
for dust temperatures of ~800 K, gas densities are only 
Ug ^ 10^ cm^'^. These densities are not prohibitively 
high and are, in fact, close to what would be expected 
for typical progenitor mass loss rates {^ 10~^ — 10~^ 
Mq yr"^) and wind velocities (^ 10 km s^^) at radii of 
- 10" cm. 

The authors also argue that coUisional heating does 
not predict the observed drop in the infrared luminos- 
ity before its X-ray counterpart. We find, however, that 
the sputtering lifetime, Tgput, of the dust grains is short 
enough that there should not necessarily be a correlation 
between the infrared and X-ray luminosities. The sput- 
teri ng lifetime a t gas temperatures above ^10^ K given 
by (|Dwek fc Ar cndt 1992) is 



' sput 



(yr) « 1 X 10^ 



a(/im) 
ig(cm^'^ 



(4) 



For the grain size (a = 0.1 fiia) and gas densities {ug ~ 
10^ cm^"^) above, Tgput ~ 1 year, which is shorter than 
the observed light plateau in SN 1995N. In such a short 
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Fig. 5. — Near-infrared fluxes for SN 2005ip at (a) early (top) and (b) late (bottom) times. Presuming a dust emission/absorption 
efficiency Q oc -^, where n is an integer, a Planck curve is flt to the H — Kg colors to establish a near-infrared dust temperature. We 
consider both large grains that are associated with condensing dust (n = 0) and smaller grains that are associated with pre-existing 
optically thin dust (n = 1 & 2). The Planck curves are plotted over the original fluxes. The curves exclude the J-band fluxes because of 
SN photosphcric flux and line emission contamination, especially at early times when the photospheric temperature is ~10,000 K. These 
fluxes consistently fall above the Planck curves. Overall, the Planck peaks trend towards longer wavelengths/cooler temperatures. 



period of time, infrared emission from coUisionally heated 
dust will decrease as the larger dust grains are destroyed. 
These effects could explain the lack of correlation be- 
tween the infrared and X-ray luminosities. None theless, 
collisional heating is still unlikely in this case. iDraind 
(|l981i) shows that coUisionally heated grains radiate only 
~ 12% of the shock power before the grains are sput- 
tered, thereby limiting the contribution of infrared emis- 
sion from coUisionally heated dust. 

SN 2006JC also shows an extended period of increased 
late-time infrared emission. Unlike the previous two ex- 
amples, SN 2006JC is a Type lb event with a dense cir- 
cumstellar medium that is likely associated with a LBV- 
like outburst originati n g from the progenitor a t -730 days 
(jYamaokaet al.ll2006l : iPastorello et all |2007| ) . Between 
days 51 and 75, the supernova undergoes a brighten- 
ing at near-infrared wavelengths (lArkharov et al.l 120061 : 
iDi Carlo et al][200l . On day 75. iSmith et'airpioi l 
measured a luminosity L ~ 2.5 x lO''^ erg s~^ and tem- 
perature T r^ 1600 K. The temperature stayed constant 
for only ~1 month and then cooled. Throughout this 
time period, X-rays indicate interaction of the supernova 
shock with the circumstellar shell cr eated by the previous 
LBV outburst (jimmler et al.ll2008D . The SN continued 
to remain bright at infrared wavelengths thro ugh day 
493, which was the last ob servation epoch (jMattila et al.l 
[200a IPTCarlo et al.l[200l . 

While warm dust is the widely accepted source of the 
late-time n ear-infrared emissio n, the origin of the dust 
is debated. ISmith et al.l (|2008l ) argue the dust is newly 
formed based on (1) the development of a red/IR con- 
tinuum consistent with the temperature of newly formed 
graphite grains or slightly hotter silicate grains, (2) the 
coincidental fading of the redshifted wings of the narrow 
He I emission lines, a nd (3) a simultane ous decline in the 
expected optical f lux (iFolev et al.ll2007D a nd brightening, 
in nea r-infrared. ISakon et al.l ()2007[ ) and iNozawa et al.l 
()2008l ) propose t he dust conde nses in the freely expand- 
ing ej ecta. Both ISmith et all (J2008) and iMattila et all 
(|2008[ ) reject this scenario because it cannot explain the 
observed spectral line reddening. Instead, these authors 
favor dust condensation in the cool, dense shell created 



as the su pernova shock inte racts with the circumstellar 
medium. ISmith et al.l (|2008f ) favor dust condensation in 
the reverse shock because of the observed carbon-rich fea- 
tures that are consistent with enric hed ejec ta encounter- 
ing the reverse shock. In contrast, IMattila et al.l ()2008l ) 
conclude the reverse shock would not be radiative and, 
therefore, not conducive to dust formation. Instead, the 
authors argue for dust condensation in the cool, dense 
shell behind the forward shock. 

In any case, the temperature remains constant for 
about one month at early times. ISmith et al.l 



gue that if there is a sufficient amount of circumstellar 
material, then grains will continue to form at the vapor- 
ization temperature. To produce the observed luminosity 
of L = (2-3) X 10^1 erg s^^ by cooling through -1600 K, 
the particle condensation rate is given as A^ = L/e, where 
e is the thermal energy per particle. For dust grains, 
the thermal energy per particle is given as e = CAT , 
where C is the dust heat capacity (Drainc fc Lil l2001h 
and AT is 1600 if we assume the dust cools immediately 
and contributes only to the instantaneous luminosity. In 
this case, we calculate that the dust must condense on 
the order of —1000 M© day""'^. This unrealistic mass 
condensation rate indicates that there must be an ad- 
ditional heating source to power the observed infrared 
luminosity. Given the optical light cu rve, which is pow- 
ered by radioactivity, auickl v dechnes (iFolev et 311120071) 
at about the same time as 'Arkhar ov et al.l ( 2006f ) ob- 
serve near-infrared brightening suggests that the newly 
formed dust is absorbing the optical emission from the 
central supernova. In this case, the likely dust heating 
source is the photo spheric emission from radioactivity. 
IMattila et al.l ()2008l ) similarly argue that the cooling of 
the newly formed grains is insufficient to account for the 
observed near-infrared luminosity and conclude a similar 
heating source for their dust condensation scenario. 

3.3. The Source of the Infrared Excess in SN 2005ip 

Warm dust surrounding SN 2005ip is the likely source 
of the near- infrared excess. As discussed in iJU there are 
several possible explanations for the large late-time Kg- 
band luminosity: dust condensation in the ejecta, col- 
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Fig. 6. — The (a) temperature (top) and (b) luminosity (bottom) evolution for SN 2005ip for all dust emission/absorption efficiencies. 
At early times (<150 days), the derived temperature drops as contamination from the SN photosphere fades. At later times (>150 days) 
the dust emission dominates. Although the exact point at which the dust temperature begins to dominate is unknown, we presume the 
photospheric component to decline sufficiently below the dust temperature after ~150 days. While the early dust temperatures vary 
with emission/absorption efficiency by a couple hundred Kelvin, they are all less than or equal to the dust vaporization temperature 
(~ 1400 — 1600 K). These near-infrared observations probe only the warmest and most luminous grains. Cooler dust may also be present, 
but its emission peaks at longer wavelengths. Assuming a distance of 30 Mpc, integrating the Planck curves over all wavelengths yields an 
approximate bolometric dust luminosity. In contrast to the decreasing dust temperature, the dust luminosity plateaus from early times. 
It is unlikely that radioactive decay within the ejecta (solid line) is the energy source as the decay curve does not match the observed 
luminosities. 



lisional heating in circumstellar medium shock interac- 
tion, an infrared echo, and dust condensation in the cool, 
dense shell behind either the forward or reverse shocks. 
This section briefly explores each possibility in order to 
distinguish between the different scenarios. SNe f998S, 
1995N, and 2006jc provide important points of reference. 

3.3.1. Dust Condensation in the Unshocked Ejecta 

Dust condensation may occur in the SN ejecta as it 
expands and cools. The available power in the expand- 
ing gas, primarily from radioactive decay, dictates the 
temperature and luminosity evolution of the grains. The 
decay of ^^Co, with a half-life of 88 days, dominates the 
plausible radioactive power sources over the first hun- 
dreds of days. Other radioactive sourc es with longer 
half- li ves become dominant at later times. IWooslev et alJ 
()1989| ) describe the details of radioactive energetics in su- 
pernovae. For SN 2005ip, the radioactive luminosity is 
larger than the observed dust luminosities at early times, 
but later falls well below the observed dust luminosity 
(Figure [6]), for an assumed 0.1 Mq of ^^Ni formed in the 
explosion. Furthermore, the observed dust luminosity 
varies slowly and does not match the radioactive decay 
rate of ^^Co. While dust condensation may have oc- 
curred in the freely expanding ejecta of SN 2005ip, the 
available radioactive power cannot produce the observed 
dust grain luminosities at late times. 

3.3.2. Collisional Heating 

As the forward shock expands into the circumstellar 
medium, the hot post-shocked electrons may collisionally 
heat pre-existing dust grains produced by the progenitor 
wind. To heat the dust to t he observed t empe ratures 
(T « 1400-1600 K), figure 1 of iPwek et all (pOM ) shows 
that the required gas densities are high, (n > 10^ cm~^), 
assuming post-shock gas temperatures on order of 10^ K. 
The sputtering effects for such high gas densities would 
significantly limit the contribution of infrared emission 
from collisionally heated dust throughout the duration 



of the observations. 

A stronger argument against this scenario is the in- 
consistency in the time lag between the SN and the ob- 
served infrared emission. The SN peak luminosity va- 
porizes all dust grains within a radius r^ that, for a dust 
emission/absorption efficiency Q\ oc A~^, i s proportion al 
to the supernova's peak luminosity Lp^ak (|Dweklll983l ). 

. ^ X 0.5 



T5 



(5) 



For a SN with a peak UV-optical bolometric luminos- 
ity of 1 X 10^" Lq with an exponential decline rate 
time-s c ale of 25 days (L* = 1.0 x ioWg-t/25d ^^^^ 
iDwekl ()1983l I1985D shows r„ « 6 x lO^^ cm (0.06 Hght 
years) for carbon-rich grains {Tevap = 1900 K) and 
3 X 10"'^^ cm (0.31 light years) for oxygen-rich grains 
[Tevap = 1500 K). While the very early evolution of SN 
2005ip is not well known, an early spectrum of SN 2005ip 
(Mo djaz, private communicatio n) and i?-band photome- 
try (jSmith et al.l 12008. in prepi ) indicate an observed lu- 
minosity of ^ 0.2 X 10^° Lq and an exponential time-scale 
of '^35 days. Considering the supernova was not detected 
until a few weeks past maximum, the actual peak is most 
likely greater by a factor of a few, which i s comparable 
to the supernova peak luminosity used by iDwekl ()1983D 
as described above. 

The size of this dust cavity implies a time lag between 
the initial explosion and the onset of infrared emission 
equal to the shock propagation time across the dust free 
radius. Using the minimum of the Ha absorption line, 
iModiaz et al.l (|2005l measure the SN expansion velocity 
to be 1.5 X 10^ cm s^^ (0.05c) at an early age. This speed 
refers to the photosphere and there is presumably higher 
velocity gas, but there is also likely to be deceleration 
resulting from the interaction. At this speed the shock 
must travel for ~450 days before encountering the dust, 
assuming the dust cavity radius for carbon-rich grains 
given above, and ^2300 days for the oxygen-rich grains. 
FigurelH however, shows a considerably shorter time de- 
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lay before a detectable infrared excess. 



3.3.3. IR Echo 



The IR echo model given by iDwekl ()1983D is tempt- 
ing to explain SN 2005ip because it predicts an extended 
luminosity plateau, as observed in Figure [Sb. Again, 
the fading photosphere is responsible for the large tem- 
peratures and rapidly evolving luminosity over the first 
150 days in Figure [6)3. Once the photosphere fades suffi- 
ciently, the dust temperature is ~1400-1600 K. This tem- 
perature is approximately the dust vaporization temper- 
ature, which is consistent with the dust lying just beyond 
the cavity predicted by the IR echo model. 

Assuming the IR echo scenario to be correct, the du- 
ration of the Ks-ha,nd plateau (^900 days) indicates 
the dust cavity radius must be r^ > 1 light year (~ 
10^® cm). Given the relationship in equation ([5]), this 
radius would imply a peak luminosity L ~ 1 x 10^^ Lq. 
This luminosity is ~ 10 times larger than the estimate 
made in § 13.3.21 Assumin g a fairly exponential decay 
(jSmith et al.l 12008. in preiJ ). the total radiated energy 
from the supernova over the first few weeks is approx- 
imately E ^ 10^^ ergs. 

Given the infrared luminosity plateau in Figure[5}D, the 
total radiated infrared energy through day 935 is approx- 
imately £;/« = 3x10'*! erg s-^x (8.1x10^ s) «2.4xl0'*9 
ergs. This energy is only a lower limit because little de- 
cline of the infrared flux has been observed. Further- 
more, this luminosity does not account for cooler dust 
contributions. With the simplified assumption that these 
grains are emitting as blackbodies, the optical depth is 
approximated as r « e+'e '^ 0.02, which, like the in- 
frared energy, is also a lower limit. This optical depth 
is consistent with the ultraviolet/optical observations of 
the supernova, which do not indicate a large amount of 
reddening {B — V — 0.2 on day 461) (jimmler fc Pooleyl 
[2P07). 

The optical depth gives a good estimate of the gas mass 
loss from the progenitor. The above optical depth r ~ 
0.02 corresponds to a visual extinction Ay = 1.086 x r = 
0.02. Assuming dust properties comparable to the Milky 
Way's interstellar medium, the column density is given as 
Nh = 1.8x10'^^ xAv =3.9xl0i^cm"2^ Given the inner 
dust cavity radius derived above, r^ ~ 10^^ cm, the total 
mass loss is then Mioss ~ N^ x {'inrl) x niH ~ 0.5 Mq. 
Assuming a wind velocity of '^10 km s~^, the mass loss 
rate is ~ 2 x 10~^ Mq yr~^, which agrees with measured 
values of typical SN winds. 

Despite the consistency in the above arguments, the 
required luminosity for this scenario is remarkably high 
for a SN and is comparable t o the most luminous events 
that have ever been observed (jQuimbv et al.l[2007h . Even 
if an eruption prior to the core collapse could have pro- 
duced the large dust cavity, the same peak luminosity 
is required to warm the dust to the observed temper- 
atures. Early time optical observations do not suggest 
that the SN was unusually brig ht at maximum (Teamo 
2007, Smith et al.'2 008. in pred ). An IR echo, while able 
to adequately predict an extended luminosity plateau, is 
not the likely source of the late-time infrared emission in 
this case. 

3.3.4. Dust Condensation in Post-Shocked Gas 



A plausible source of warm dust in SN 2005ip is dust 
condensation in the shocked region that results from the 
collision of the ejecta with the pre-existing circumstellar 
medium. If the shocked gas is able to efficiently cool, 
the high pressure in the region results in gas densities 
> 10^° cm~^ when the temperature drops to ~ lO'^ K. 
The shocked region is bounded by a forward shock in the 
circumstellar gas and a reverse shock in the supernova 
ejecta. Radiative cooling is more likely at the reverse 
shock because of the higher density, lower shock velocity, 
and the possibility of heavy element enrichment. Cooling 
at the forward shock is less likely, but the presence of 
clumps in the circumstellar medium, as indicated by the 
relatively narrow lines observed in SNe Iln, may allow 
radiative shocks. 

While we cannot provide any direct spectroscopic evi- 
dence at this time to support dust condensation in a par- 
ticular shock front, the available photometry is consistent 
with dust condensation in the shocked region. The ob- 
served dust temperatures and luminosities (both infrared 
and X-ray) are similar to SNe 1998S, 1995N, and 2006jc. 
Like S N 2005ip, these three supernoyae have associated 
X-ray s (|Fox et al.ll2000l : lPoolev et al.ll2002t [immler et all 
l2008f ). which are evidence for interaction of the supernova 
shock with the circumstellar gas. These results indicate 
that the source of the infrared emission is most likely 
associated with the shock interaction region. While SN 
1995N was the only one of these supernovae not previ- 
ously associated with dust condensation in the reverse 
shock, we consider the possibility below in ^ 

The dust temperature (T « 1400 - 1600 K) on day - 
200 is also consistent with the condensation temperature 
of carbon-rich grains. Assuming the dust shell emits as 
a single, spherical blackbody (r > 1), a minimum radius 
of ~ 10^^ cm is consistent with the observed luminosity. 
This radius corresponds to a shock expansion velocity of 
^ 3000 — 8000 km s~^, where the range is related to the 
grain composition. The upper part of this range is quite 
reasonable. Unlike the IR echo scenario, this picture is 
physically coherent and does not require an unreasonably 
high luminosity. 

The dust equilibrium temperature is ultimately de- 
termined by a balance of heating and cooling. In this 
case, the dust heating mechanism is likely radiative shock 
emission. While the observ ed X-ray lum inosity on day 
466 is only 1.6 x 10^° erg s~^ (jimmler fc Po olcv 2007), the 
corresponding optical luminosity is ~ 2 x 10^"'^ erg s~^, 
which is deduced from the V magnitude at this time. 
This luminosity is comparable to the observed infrared 
luminosity (r « r^^'J^ ~ 0.5). The i?-band luminos- 

■' } E+Ejft J. ^ 

ity plateau ([Smith et al.ll2008. in prepD suggests that this 
visible component persists throughout the duration of 
our observations and, if absorbed by the newly formed 
dust, is a significant source of heat for the dust grains. 

As the cool, dense shell expands, the radiation field 
from the shock emission becomes diluted and we ex- 
pect the dust temperature to vary as T ex (L/r^)^/'*. 
Assuming an expansion velocity of r^ 5000 km s~^, by 
day ^ 800 the dust shell radius would have expanded 
by about a factor of 3 and the temperature would have 
dropped from ~ 1500 K to ~ 850 K. Figure [Sh shows 
a decreasing T with time, although the decline appears 
to be slightly weaker than expected. This discrepancy is 
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most likely related to the accuracy of our assumptions 
involving expansion velocity and grain composition. 

4. DISCUSSION AND CONCLUSIONS 

This paper presented near-infrared observations of SN 
2005ip for the first 935 days following detection. A large 
Ks-hand luminosity persists even as the supernova's J- 
band luminosity falls. Among a variety of potential 
mechanisms, dust condensation in the cool, dense shell 
downstream from the reverse shock is the likely source of 
the observed infrared emission. 

We are able to rule out other mechanisms for late- 
time near-infrared emission. While dust condensation 
in the supernova ejecta is a likely possibility, the du- 
ration of the observed near-infrared light curve is in- 
consistent with the quickly declining radioactive heating 
source. Shock/mechanical heating of pre-existing dust 
grains from prior mass loss is not possible because the 
near infrared excess appears quickly (<100 days), provid- 
ing insufficient time for the expanding ejecta to cross the 
dust-free cavity formed by the supernova. An IR echo 
successfully explains the fairly uniform luminosity over 
the period of observation. The duration of the infrared 
excess, however, implies a cavity size that requires a peak 
supernova luminosity that is much larger than observed 
for SN 2005ip and is comparable to the most luminous 
supernova ever observed. 

On the other hand, there is mounting evidence to sup- 
port dust condensation in the shock interaction region. 
SN 2005ip shares many properties with SNe 1998S and 
2006JC, for which spectra reveal dire ct evidence of dust 
condensation in the reve rse shock (jPozzo et al.l l2004t 
ISmith et al.ll2008. in prenf ). While we currently have no 
such spectra for SN 2005ip, the observed infrared lumi- 
nosities and temperatures are consistent with dust form- 
ing in shocked, cooled gas. The observed X-rays indi- 
cate the infrared emission is most likely associated with 
the shock interaction region. Furthermore, the initial 
dust temperature of ~ 1400 — 1600 K is consistent with 
the dust condensation temperature of carbon-rich grains, 
which we expect to find in the ejecta encountering the 
reverse shock. The temperature evolution is also some- 
what consistent with the decline predicted by an expand- 
ing cool, dense shell and diluted radiation field. Unlike 
the IR echo, this scenario is physically plausible and a 
large initial luminosity is not required. 

Of course, near-infrared photometry is only a small 
part of a much larger picture. Infrared spectra can 
better constrain the dust temperature, infrared flux, 
and emission/ absorption efficiency than can be obtained 
with near-infrared photometry alone. Mid-infrared spec- 
tra also hold the prospects of revealing the mineralogy 
of the heated dust, distinguishing between carbon-rich 
(graphite) and oxygen-rich (silicate) grains. Based on 
the dust composition and mass loss history, the nature 
of the progenitor system may be deduced. Optical spec- 
tra are also of interest in that they can potentially reveal 
more evidence for dust fo rmation and heating via cir- 
cumstellar interaction (e.g. iFransson et al.l[200a ) 

The overall similarity between SN 1995N and SN 
2005ip suggests that SN 1995N be considered in the 
context of the dust condensati on model, as o ppose d 
to the echo model favo r ed by iGerardv et al.l (|2002D . 
Indeed, iPastorello et al.l (|2005[ ) show that the late- 



time optical and infrared luminosities of SN 1995N are 
comparable, which suggests a correlation between the 
ejecta/circumstellar medium interaction and the warm 
dust. The IR echo model also predicts a optical light 
echo from the scattered light, as well as an infrared echo 
produced by the absorpti on of the supernova peak lu- 
minosity (|Chevalierlll986f ). Yet detailed spectroscopic 
observations of SN 1995N suggest the late visible emis- 
sion is due to circum stellar interaction and not an echo 
(JPra nsson et al.ll200 2l . These spectra reveal evidence for 
absorption in the diminution of the red sides of emis- 
sion lines compared to the blue sides, as would be ex- 
pected with dust formation. Furthermore, the authors 
present evidence for an increased helium abundance, 
CNO-burning products, and Lya-pumped fluorescence of 
Fe II lines, all of which suggest a complex interaction be- 
tween the ejecta and a pre-existing, dense circumstellar 
medium. The dust temperature of SN 1995N on day 730, 
700-800 K, is slightly smaller than for SN 2005ip. While 
this difference may be not be statistically significant, it is 
possible that SN 1995N had either a significantly higher 
expansion velocity or a different grain composition. 

Our discussion of SN 1995N and the data for SN 2005ip 
indicate that Type Iln SNe may be sources of dust in our 
universe. This association would explain the consistent 
late-time in frared emission observed in several Type Iln 
events (e.g. IPastorello et al.l 120021: IGerardv et all l2002t 
iPozzoet all |2004[) . The rarity of Type Iln events, 
which represent only ~2-3% of all core-collapse SNe 
(jGal-Yam et al.M2007| ). makes collecting complete light 
curves difficult. Nonetheless, an interesting study would 
include follow-up observations of all Type Iln events over 
the past several years to determine which events are still 
bright in the i nfrare d. We also reiterate the call by 
IGerardv et al.l ()2002f ) for simultaneous UV and X-Ray 
observations to better constrain r and the amount of 
shock radiation being reprocessed by the dust grains. 
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